Stability of the r-modes in white dwarf stars 
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Stability of the r-modes in rapidly rotating white dwarf stars is investigated. Improved estimates of 
the growth times of the gravitational-radiation driven instability in the r-modes of the observed DQ 
Her objects are found to be longer (probably considerably longer) than 6 x 10 9 y. This rules out the 
possibility that the r-modes in these objects are emitting gravitational radiation at levels that could 
be detectable by LISA. More generally it is shown that the r-mode instability can only be excited in 
a very small subset of very hot (T> 10 6 K), rather massive (M > O.9M0) and very rapidly rotating 
(fmin < P < 1.2P m i n ) white dwarf stars. Further, the growth times of this instability are so long 
that these conditions must persist for a very long time [t > 10 9 y) to allow the amplitude to grow to 
a dynamically significant level. This makes it extremely unlikely that the r-mode instability plays 
a significant role in any real white dwarf stars. 



I. INTRODUCTION 

Recently it was discovered by Andersson jjj , and Fried- 
man and Morsink |g| that gravitational radiation tends 
to drive unstable the r-modes of all rotating stars. Lind- 
blom, Owen, and Morsink || subsequently calculated the 
strength of the gravitational radiation coupling to these 
modes and found it to be strong enough to overcome vis- 
cous dissipation in hot rapidly rotating neutron stars. 
This calculation was confirmed using better estimates 
of the dissipative coupling by Andersson, Kokkotas, and 
Schutz U and Lindblom, Mendell, and Owen ||. Thus 
it is now generally expected that this r-mode instability 
will limit the rotation rates of hot young rapidly rotating 
neutron stars. The excess angular momentum in these 
stars will be carried away by gravitational radiation, and 
these stars will be spun down to more slowly rotating sta- 
ble configurations within about one year of their births. 
The gravitational radiation emitted during this process 
may eventually be detectable by LIGO, (see Owen, et 
al. ||, and Brady and Creighton Q). 

The possibility that this r-mode instability plays a role 
in other astrophysical systems has also been proposed. 
For example, the possibility that this instability might 
play a role in limiting the angular velocities of old and 
relatively cold neutron stars spun up by accretion was 
considered by Bildsten M, Andersson, Kokkotas, and 
Stergioulas |J] , and Levin [|l0fl . It has also been proposed 
that this instability might play a role in rapidly rotat- 
ing white dwarf stars by Andersson, Kokkotas, and Ster- 
gioulas II and by Hiscock |ll ] . This last possibility will 
be investigated more thoroughly in this paper. A set of 
rapidly rotating white dwarf models, and the timescales 
for the gravitational-radiation driven r-mode instability 
in these stars are computed in Sec. [n| The possibility 
that this instability is playing a significant role in the 
observed DQ Her objects is examined in Sec. III. It is 



shown that the minimum growth time for this instability 
in these objects is 6 x 10 9 y, and that the actual growth 
time is almost certainly much longer than this. Thus, 
the r-mode instability is not playing any significant role 



in these objects, and gravitational radiation from the r- 
modes in these objects will not be observable by LISA. 
Sec. evaluates the effects of viscosity on the r-modes 
in white dwarf stars. It is shown that all white dwarfs 
with core temperatures cooler than about 10 6 K are sta- 
ble; that all white dwarfs with masses less than 0.9M Q 
are stable; and that all white dwarfs with rotation pe- 
riods longer than 1.2P m i n (where P min is the minimum 
rotation period of that star) are stable. It is shown that 
white dwarf stars cool too quickly (in the absence of ac- 
cretion) for the r-mode instability to grow to significant 
levels in any star. And finally, it is shown that this in- 
stability can only play a significant role in white dwarf 
stars that are maintained at very high core temperatures 
by accretion, and that remain at nearly maximal rota- 
tion rates for a period of time that exceeds about 10 9 y. 
It seems very unlikely that these conditions will ever be 
met in any real white dwarf stars. 



II. RAPIDLY ROTATING WHITE DWARFS 

In order to study to properties of the r-modes in white 
dwarf stars, simple models have been constructed using 
the equation of state of a zero-temperature degenerate 
Fermi gas with Z/A = i that is appropriate for carbon- 
oxygen white dwarfs JL2|. Families of rigidly rotating 
stars were constructed using the numerical techniques 
described by Ipser and Lindblom |l3f| for finding very 
rapidly rotating and highly non-spherical models. A con- 
stant mass family of rotating models was constructed by 
successively spinning up a non-rotating model until its 
angular velocity was equal (or very nearly equal) to the 
frequency of the equatorial orbit located just at the sur- 
face of the star. In Fig. [I] are presented the rotation 
periods P m in of these maximally rotating white dwarf 
models as a function of the mass of the star. Also shown 
in Fig. [l] are a set of points that represent the minimum 
rotation periods for white dwarf stars as determined by 
Hachisu H|. The numerical models constructed here 
used a spherical grid that allowed a much finer determi- 
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FIG. 1. Minimum rotation periods (in seconds) for rigidly 
rotating white dwarf stars based on the equation of state of a 
zero-temperature degenerate Fermi gas with Z/A = i. 
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FIG. 2. Gravitational radiation growth times (in years) of 
the r-modes for white dwarf stars rotating at their maximum 
possible rotation rates. 



nation of the location of the surface of the star than the 
models constructed by Hachisu. Nevertheless, there is 
good agreement between our calculations of the limiting 
rotation periods of these stars. 

The growth time of the gravitational-radiation driven 
r-mode instability is estimated using the expressions de- 
rived by Lindblom, Owen, and Morsink ||. In particular 
the growth time t gr for the dominant m — 2 mode is 



t gr 25 \ 3j c J 1 \ c 1 



(2.1) 



where p and f£ are the density and angular velocity of the 
equilibrium stellar model, and G and c are Newton's con- 
stant and the speed of light respectively. This is only the 
lowest order term in the expansion for 1/t gr in powers of 
the angular velocity. To lowest order therefore it is suffi- 
cient to evaluate this integral using the density function 
p for the non-rotating model of a given mass. The grav- 
itational radiation growth time t gr defined in Eq. (2.1) 
is proportional to fl~ 6 . It is convenient to define a char- 
acteristic growth time f GR therefore for the star rotating 
at the maximum possible angular velocity. In particular 
then the gravitational radiation growth time is related to 
this characteristic growth time by 



T Gf 



P 



(2.2) 



The values of the characteristic growth time f GR for a 
range of white dwarf masses are presented in Fig. 



III. DQ HER OBJECTS 

Recently Andersson, Kokkotas, and Stergioulas Q and 
Hiscock |ll]] have suggested that the instability of the r- 
modes driven by gravitational radiation may be playing 



an important astrophysical role in the DQ Her type cat- 
aclysmic variable systems. These objects emit radiation 
with periodic luminosity variations that have been inter- 
preted as the rotation periods of rapidly rotating white 
dwarf stars. The systems having the shortest periods 
are WZ Sge with period 27.9s, AE Aqr with 33.1s, V533 
Her with 63.6s, and DQ Her with 71.1s U|. If these 
represent the rotation periods of the white dwarfs (a n 
interpretation that is not universally accepted []il|[l7|| ) 
then Fig. [I] shows that these systems could (depending on 
their masses) contain very rapidly rotating white dwarf 
stars. The gravitational radiation growth times of the r- 
modes in these systems are smallest — for a star of given 
mass — for the most rapidly rotating stars. It can also be 
shown (by brute force numerical examination of the data 
in Fig. |2| together with Eq. [T^) that the gravitational ra- 
diation growth times of the r-modes in these systems are 
also smallest — for stars of given rotation period — for that 
star which has the smallest mass. Hence it is straight- 
forward to show from the data in Figs, [j] and || that the 
minimum possible gravitational radiation growth times 
for the r-modcs in the DQ Her objects are 6.7 x 10 9 y for 
WZ Sge, 1.5 x 10 10 y for AE Aqr, 3.7 x 10 n y for V533 
Her, and 6.6 x 10 n y for DQ Her. 

The estimates of the gravitational radiation growth 
times given above are almost certainly significant un- 
derestimates. First, internal fluid dissipation (i.e. vis- 
cosity) has been completely neglected. This will always 
act to increase these timescales. Second, the assump- 
tion that these are maximally rotating white dwarf stars 
is quite questionable and consequently the gravitational 
radiation growth times are expected to be much longer 
than the values given above. For example, the most re- 
cent mass determination of the white dwarf in WZ Sge is 
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only 0.3M Q 0.0 The minimum gravitational radiation 
growth time for a star of mass O.3M is 1.6 x 10 12 y. The 
most recently measured mass for AE Aqr is O.8M [Q, 
and for DQ Her is 0.6M Q JnJ. These values imply that 
these objects are rotating significantly below their max- 
imum rotation rates: P/Pmin — 1-4 for AE Aqr and 
P/Pmin = 2.0 for DQ Her. The gravitational radiation 
growth times for these systems would then be 2.4 x 10 10 y 
for AE Aqr and 1.3 x 10 12 y for DQ Her. And further, re- 
cent observations indicate that the rotation period of the 
white dwarf in DQ Her is 142.2s rather than 71.1s [ |l6| . 
This increases the gravitational growth time for this star 
to 8.1 x 10 13 y. 

In order for the instability in the r-modes to play 
a significant dynamical role in the evolution of a sys- 
tem, the dimensionless amplitude of the mode must grow 
to a value of order unity. For example, the gravita- 
tional radiation amplitudes computed by Hiscock ]Tl| ] 
are based on a presumed balance between accretion and 
gravitational-radiation torques on the star. Using the 
equations in Owen et al. H, it is easy to show that (in 
a rapidly rotating star) this balance requires the dimen- 
sionless amplitude of the r-mode to be a constant of order 
unity multiplied by (t gr M /M) 1 / 2 . Given the observed 
accretion rates, the amplitudes of the r-modes would 
have to be of order unity in these systems to maintain 
this balance. The extreme lengths of the gravitational- 
radiation timescalcs in these systems means that there 
is not enough time for the amplitude of the r-mode to 
grow to such a large value. So even if viscosity were 
unimportant and all white dwarfs were unstable to the 
gravitational-radiation driven r-mode instability, this in- 
stability can not be playing any significant role in the 
observed DQ Her objects. Gravitational radiation from 
the r-mode instability in these objects will not be de- 
tectable by LISA. 



where tv is the viscous timescale. For the m — 2 r- 
mode of primary interest to us here, r v is given by the 
expression pj 




r/r^dr 



(4.2) 



where rj is the shear viscosity. The dominant form of 
shear viscosity in hot white dwarf stars is expected to be 
electron scattering with the ion liquid. Nandkumar and 
Pethick [^0| provide the following analytical fit for the 
shear viscosity in a 12 C ion liquid, 
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where I2 is 
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The density pe and temperature T§ are to be given here in 
units of 10 6 gm/cm 3 and 10 6 K respectively. This formula 
for 77 is sufficiently accurate (i.e. to within about 10%) 
for the densities above 10 6 gm/cm 3 wh ich dominate the 
integral in the denominator of Eq. ( fLlq) 



The growth time r defined in Eq. (4.1) is negative for 
slowly rotating stars: P ^> P m in- In this case viscosity 
dominates and the r-mode instability is suppressed. For 
more rapidly rotating stars, however, this expression may 
become positive and hence the r-modes unstable. The 
critical rotation period P c where the instability first sets 
in is determined by setting 1/r = 0. It follows then that 



IV. VISCOUS DISSIPATION 

Shear viscosity tends to suppress the gravitational ra- 
diation driven instability in the r-modes. In particular 
the growth time r of the mode becomes 
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(4.1) 



*A white dwarf with rotation period 27.9s must have at min- 
imum a mass of 0.725Mq. Therefore, if the rotation period of 
WZ Sge is 27.9s, then its mass must be greater than O.725M0 . 
Conversely, if its mass is less than O.725M0 as suggested by 
recent mass determinations, then the observed 27.9s period 
must be a harmonic of its fundamental rotation period, or be 
caused by a pulsation or other phenomenon unrelated to its 
rotation. 
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This critical rotation period P c has been evaluated for the 
white dwarf models discussed in Sec. [0] and a range of 
core temperatures appropriate for hot white-dwarf stars. 
These results are depicted graphically in Fig. These 
results show that no white dwarf is unstable to the r- 
mode instability when its core temperature falls below 
10 6 K. Further, no white dwarf star with mass smaller 
than 0.9Mq is subject to the r-mode instability if its 
core temperature is below 2 x 10 8 K, the maximum core 
temperature expected to occur in accreting white dwarf 
systems ||l],^]. No white dwarf with rotation period 
greater than 1.2P m j n is unstable to the r-mode instability 
if its core temperature is below 2 x 10 8 K. 

To investigate further whether the r-mode instability 
might be playing a role in some white dwarf stars, the 
growth times r for these modes are examined. The ex- 
pression for this timescale, Eq. (4.1), can be re- written 
in terms of the quantities graphed in Figs. |l|, || and [| 
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FIG. 3. Critical rotation periods P c where the gravitational 
radiation instability in the r-modes first sets in: stars with 
P < P c are unstable. The mass of the white dwarf (in units 
of Mq) labels each curve. 
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FIG. 4. Minimum growth times r m i n for the gravitational 
radiation driven instability of the r-modes for white dwarf 
stars. The mass of the star (in units of Mq) labels each 
curve. The dashed curves depict the time r c needed for a 
white dwarf to cool to a given core temperature. 
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The timescale r becomes infinite as P — > P c and is mini- 
mum for P = P m i n . Thus it is useful to define the small- 
est possible growth time r min for a star of given mass and 
temperature: 
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The values of these minimum growth times are depicted 
in Fig. [I| It follows that no white dwarf has an r-mode 
instability with growth time shorter than 7 x 10 7 years if 
its core temperature is below 2 x 10 8 K. 

In the absence of accretion a white dwarf star will 
quickly cool. Its core temperature drops to the value T 
in approximately the time r c given by the expression [^3| , 
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The dashed curves in Fig. [| represent this cooling time 
t c and also 0.1t c . It follows that most of the minimum 
r-mode instability growth times are longer than r c . In 
these stars, the amplitude of the r-mode will grow by 
less than a factor of e before the star cools and the insta- 
bility is suppressed. No star has a minimum instability 
growth time that is shorter than 0.1r c . Thus, the r-mode 
instability will never have time to grow to a dynamically 
important level in any white dwarf except during the time 
it is heated by accretion. 

If a white dwarf is heated by accretion, this heating 
can last for only a finite time, given approximately by 
t a = AM/M where AM is the total amount of accreted 
material and M is the accretion rate. Since the r-mode 
instability is never effective for stars with M < 0.9M Q 



we see that the r-mode instability can only act for a 
time when AM < 0.6M©. Further, we see from Fig. || 
that T min > 7 x 10 7 y for stars with T < 2 x 10 8 K, 
the maximum core temperature reached in models of 
accretion onto white dwarfs pTjJ?2| |. In order for the 
r-mode instability to have time to grow to a dynami- 
cally significant level then, it would be necessary to have 
t a > 10r m i n . This implies that the accretion rate must 
satisfy M < 9 X 10~ 10 MQy~ 1 . Given this low accretion 
rate, this argument can be strengthened slightly. For 
white dwarfs with M < 10" 9 M Q y- 1 the core temper- 
ature is only expected to reach about 8 x 10 7 K [gl 24 
during accretion and so from Fig. [| r m i n > 10 8 y. For 
stars with core temperatures in this range, Fig. || im- 
plies that instability can only occur on a timescale fast 
enough to act effectively within the age of the universe 
for stars in the mass range 1.05 < M < 1.4M Q and so 
AM < 0.35M Q . Thus, the upper limit on the accretion 
rate can be improved slightly to M < 3.5 x lO _lo M0y _1 . 

One further difficulty to achieving a dynamically sig- 
nificant r-mode instability in white dwarf stars is illus- 
trated in Fig. |5|. The minimum growth times T m i n shown 
in Fig. ^ only pertain to stars that are maximally rotat- 
ing P = P m in ■ The actual growth times of these modes 
increase rapidly for more slowly rotating stars, as illus- 
trated in Fig. |]for the 1.3M Q stellar model. For example, 
the growth time increases from about 10 8 y for the maxi- 
mally rotating model with T = 2 x 10 8 K to about 10 9 y for 
the still very rapidly rotating model with P = 1.15P m i n . 
Thus the very hot (T > 3 x 10 6 K) and very rapid-rotation 
conditions (P m in < P < l-15P m m) needed to allow the 
growth of the r-mode, must be maintained for a very long 
time. This can only be achieved by a very low and steady 
accretion rate that adds (or removes) angular momentum 
from the star at just the rate needed to maintain max- 
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FIG. 5. Rotation periods P (as a function of temperature) 
for 1.3Mq white dwarf stars where the growth time of the 
m — 2 r-mode has the prescribed value. The dashed curve 
is the critical rotation period P c where the growth time is 
infinite. 

imal rotation as the star's mass increases.^ In the most 
favorable case (i.e. by keeping temperatures above 10 8 K, 
rotation periods less than 1.05P m i n and the mass of the 
star near 1.3M ) these special conditions need last only 
about 10 9 y in order to allow the amplitude of the mode 
to grow by a factor of e 10 «2x 10 4 . In the more typical 
and less favorable cases this type of accretion would be 
needed for 10 10 y or longer to achieve this amplification 
of the mode. It seems rather unlikely that these special 
conditions have ever been met in any real white dwarf 
stars. 
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